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We study the behavior of an assembly of N granular particles contained in two compartments 
within a simple kinetic approach. The particles belonging to each compartment collide inelastically 
with each other and are driven by a stochastic heat bath. In addition, the fastest particles can 
change compartment at a rate which depends on their kinetic energy. Via a Boltzmann velocity 
distribution approach we first study the dynamics of the model in terms of a coupled set of equations 
for the populations in the containers and their granular temperatures and find a crossover from a 
symmetric high-temperature phase to an asymmetric low temperature phase. Finally, in order 
to include statistical fluctuations, we solve the model within the Direct Simulation Monte Carlo 
approach. Comparison with previous studies are presented. 

PACS numbers: 02.50.Ey, 05.20.Dd, 81.05.Rm 



INTRODUCTION 



In recent years considerable progress has been achieved in the understanding of granular gases, i.e. assemblies of 
moving inelastic solid grains mutually colliding and loosing a little energy in each collision 0, 0] ■ They exhibit a 
fascinating and rich phenomenology which comprehends clustering 0, 0, |l| spontaneous vortex formation Q , and 
breakdown of kinetic energy equipartition in granular mixtures 0, la] ■ More fundamentally granular gases represent 
one of the prototypical non equilibrium systems. 

Major progress have been achieved by considering the simplest situations, i.e. spatially homogeneous and time- 
independent systems. In spite of that, inhomogeneities do not always represent a nuisance or an undesired complica- 
tion, but, on the contrary, they can be a source of new insight. Within such a perspective, external fields have been 
introduced on purpose in order to probe the inhomogeneous behavior of granular gases. The seminal work of Schlicht- 
ing and Nordmeier |9( has stimulated a vivid interest in the behavior of the so called compartmentalized systems. 
A box, whose base moves periodically up and down, is separated by a vertical barrier into two compartments which 
communicate through an orifice. For strong shaking the two halves are equally populated, as it would occur in the 
case of a standard molecular fluid, whereas for weaker driving the symmetry is spontaneously broken. The resulting 
scenario resembles that of an equilibrium thermodynamical phase transition, with the population difference between 
the two compartments playing the role of the order parameter and the driving intensity that of temperature. Our 
understanding of the problem has increased since then due to a series of studies. These range from new experiments 
UH El j phenomenological flux models [ToL FLU , hydrodynamic equations and Molecular Dynamics simulation . 

The scope of the present contribution is to show that by choosing a simplified yet significant model of compart- 
mentalized granular gas it is possible not only to derive a set of equations describing the evolution of the macro-state 
of the system, but also to obtain information about its microscopic fluctuations. The present approach represents a 
bridge between the phenomenological level of refs. 0,E3,0] and the statistical mechanical level [l5|. 

The present paper is organized as follows: in the first section we define the model and introduce the statistical 
description of the system, based on a Boltzmann equation for the distribution functions, modified to take into account 
the stochastic driving. At this stage we follow the strategy of integrating out microscopic fluctuations going from a 
microscopic description based on the Boltzmann equation to a macroscopic level, where only the occupation numbers 
and the granular temperatures of the two compartments, i.e. the first two moments of the distribution function are 
taken into account. This is equivalent to neglect inhomogeneities of the system at scales smaller than the linear size 
of the compartments. To obtain qualitative insight, we first try a Gaussian solution. In other words, we assume 
that the velocity PDF's are Gaussians, whose normalizations and variances (related to the occupation numbers of 
the compartments and to their granular temperatures, respectively), can be determined by means of a set of self- 
consistent differential equations. The analysis shows the existence of a transition from a symmetric phase, where 
the compartments are equally populated and are at the same granular temperature, to an asymmetric phase, where 
the two compartments have different properties. Various predictions are made: we locate the bifurcation point, 
obtain relations between the asymptotic values of the temperatures and densities in the compartments, estimate the 
characteristic times to observe the symmetry breaking, etc. In the second section, we relax the Gaussian hypothesis 
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about the nature of the velocity fluctuations and let the populations in the two compartments to fluctuate as well. 
This is done by solving numerically the full Boltzmann equation by means of the Direct Simulation Monte Carlo 
(DSMC) 0. We notice deviations with respect to the treatment of section I, in particular near the critical point. 
Finally in the last section we present our conclusions. 



Let us consider a system, composed of two compartments A and B of the same volume, Va = V/2, and containing 
Na and Nb particles, respectively. Particle pairs belonging to the same compartment may collide inelastically loosing 
a fraction of their kinetic energy, but conserving their total momentum. The post-collisional velocities (v^ , ) are 
determined by the transformation: 



where V12 = vi — V2 and a is the restitution coefficient. 

In addition, we assume that all the grains are subject to the action of an external stochastic driving force and their 
motion between two successive collisions is described by the following Ornstein-Uhlenbeck process: 



where Tf, is a measure of the intensity of the driving. Notice that in the elastic case (a = 1), the average kinetic 
energy per particle moving in a d-dimcnsional space is K = dTb/2, which is just the ideal gas value. The simultaneous 
presence of frictional dissipation and random "kicks" renders the kinetic energy of the system stationary even in the 
absence of collisional dissipation. 

To complete the model, we allow the particles contained in A (B) to move into B (A) with a probability per unit 
time, t s -1 , provided their kinetic energy exceeds a given threshold, let say T s = jjnu 2 s . Such a mechanism schematizes 
the jump process by which particles may pass from one compartment to the other by overcoming a vertical barrier of 
height h only if ^mv 2 > nigh. 

We shall make the key assumption, dictated by mathematical convenience, of ignoring the spatial gradients which 
characterize experimental situations. To be specific we stipulate that within the spatial domain, that represents a 
compartment, the system is homogeneous. The effect of the dividing vertical wall is schematized by a selection rule 
of probabilistic nature which allows only the more energetic particles to cross the barrier. 

To make analytic progress we shall transform the stochastic evolution equations for the velocities of the particles into 
a system of deterministic equations for the distribution functions. In order to achieve this description, we shall treat 
the collisions at the level of the Boltzmann molecular chaos approximation. This approximation, widely employed 
even in the case of granular systems, is equivalent to neglect correlations among the colliding particles. 

In order to study the statistical evolution of the system we assume that the single particle phase-space distribution 
function /(r, v, t) satisfies the following properties: 



THE MODEL 




Vl — |(1 + a)(vi2 ■ <x)<r 



(1) 
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where NA(B){t) is the average number of particles in compartment A (B) at instant t. The average kinetic energy per 
particle, the granular temperature Ta(b)i IS defined as: 




(5) 



In the following we shall make the assumption that to describe the essential properties of the system, we do not 
need the detailed information contained in /(r,v,i), but these can can be captured by the following coarse grained 



distributions obtained by eliminating the r dependence of the original distributions 



V A (B) Jv a 

The change of /a(b)(v, t) is given by: 



f A (B)(y,t) = ^— [ dr/(r,v,<) (6) 

VA(B) JV 4 ,m 



dtf A (yi,t) = HJaJa) + — U- /A(vi,t) + _— / A ( Vl ,t) - -0(|v x | - u,)[/ A (vi,t) - / B ( Vl ,t)] (7) 

t& \ dvi / "n, dvi r s 

where 0(x) is the Heaviside function. The first term represents the change of the distribution due to collisions, the 
second and the third terms are due to the interaction with the heat-bath and the last describes the population change 
due to the particles migrating from one compartment to the other 16] . 

In order to obtain an explicit expression for the collision integral /(/, /) in the case of inelastic hard spheres, we 
shall follow closely the derivation of Ernst and Van Noije and set: 

KJaJa) =<J d - 1 Jdw 2 J d*(v ia .a-)|^/ A (vr,t)/ A (vr,t)-/A(v 1 ,t)/A(va,t)| (8) 

The prime on the <x integration enforces the condition v 12 • <x > 0, where <r is a unit vector along the line of centers 
of the colliding spheres at contact, whereas v** represent the prccollisional velocities, which can be computed by 
inversion of eq. J2J- 

By integrating over the spatial coordinate and over the velocity it is straightforward to obtain the equation for the 
rate of change of the occupation numbers Na and N B - It reads: 



dN A {t) 
dt 



/"dvil/^vx^-Mvi^Mlvrl-^). (9) 



We see that when the typical time r s diverges, there is no particle exchange between the two compartments, thus 
Na,Nb are constant, and the temperatures reach a constant value, which depends on the heat bath properties, the 
physical properties of the particles and on their number. 

In order to obtain the energy equation, instead, we multiply both sides by rrv\r\ and integrate over coordinate and 
velocity space: 

^M&A = * Nm _ TA(t)] + WM«)T A (t) (10) 

dr T& a Va 



Va 
T s d 



dvimv?[/ A ( Vl ,t) - f B (yi,t)]0(\vi\ - u s ) (11) 
where we have introduced, following ref. [l7|. the non dimensional quantity /i 2 : 



J &v 1 vll{f Al fA)= (12) 



'^iV 2 



J dviv 2 J dv 2 J do-(vi2 ■ 6-) |-^/ A (v*" 



with Va(b) = 2TA( B )/ m - ln the following we shall set m = 1. A more compact notation is achieved if one introduces 
the collision frequency uiA{t) and the non dimensional spontaneous cooling rate 7: 

N A d _i 

UJ A = —f=—<J v A (13) 

V27T VA 
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with fid = Y{d/2) ■ ^ ne sur f ace area °f a d-dimensional unit sphere. In the remainder of the paper we shall measure 
the temperature in units of T s , time in units of t& and length in units of a. 

To proceed further analytically we choose d — 2 and assume an approximate Gaussian form of the distribution 
function: 

N A (B) 1 v 2 
f A(B) (v,t) = tt^ 1 ^ exp(-— ) (15) 

In this case one can evaluate explicitly the collision term [17| : 

>* = {i - a2) 7W(k (16) 



We find: 



7aM = (T (l-cO-^A/— (17) 
2V A V m 

= I[JV Be - T -/^ - N A e^ A ] (18) 

where the temperature T s is given by T s — ^mu 2 s . Notice that the right hand side of the equation above represents the 
difference between the incoming flux and the outgoing flux of compartment A. A similar expression has been chosen 
by Droz and Lipowski 14] on phcnomcnological grounds. The model we study allows us to obtain self-consistcntly 
the temperatures T A and T B , a feature which was not present in most of the studies dedicated to compartmentalized 
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systems |l9j. In fact, we can write: 

d[N A (t)T A (t)] = _1 [Na{Ta + Ts)e -T s /r A _ Nb{Tb + Ts)e -T s /r B] _ 2lLUA N A T A + —N A (Tb - T A ) (19) 
at t s n 

d[N B (t)T B (t)} = „1 [Nb(Tb + Ts)e -T s /T B _ Na{Ta + Ts)e -T s /T A] _ 2iuJbNbTb + l Ng{Tb _ Tb) (20) 
at T s n 

Notice that the Arrhenius type of behavior of the transition rate T s e T ^ TA results from the combination of two 
assumptions: a) the fact that only particles whose energy exceeds the threshold T s can overcome the barrier; b) the 
Gaussian ansatz for the velocity distribution functions. The latter ingredient, although very convenient for numerical 
work, could result too crude in some physically relevant situations, because the dynamics represent a severe probe 
of the extreme value statistics of the system. It is well known that the velocity distribution of driven and undriven 
granular assemblies might be characterized by fat velocity tails [16(. 

MEAN FIELD ANALYSIS 

It is convenient to rewrite the coupled equations as: 

= L [NBe -Ts/T B _ NAe -T s /T A] (21a) 

at t s 

N A ^j^ = — — [2 [N A T A e~ T "l TA - N B T B e- T °/ T *) + (N A e~ T '/ TA - N B e- T °/ T »)(2T s - T A )] 

(21b) 

- 2^ A N A T A + -N A (T b - T A ) 

n 



N dTB{t) 

B ~dT~ 



: [2{N B T B e- T '/ TB - N A T A e~ T °l TA ) + {N B e- T °' TB - N A e~ T '' TA ){2T s - T B )] 

2 (2lG) 
2 1 u B N B T B + — N B (Tb - T B ) 

n 
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FIG. 1: Variation of the eigenvalues of the dynamical matrix with respect to the granular temperature. Notice that below 
the critical temperature the less negative eigenvalue displays a maximum. This is determined by the competition between the 
inhibiting factor due to the wall and the collisional dissipation which favors a symmetry breaking. 

Let us observe that in the temperature equations there appear two types of fluxes: the first due to the unbalance 
of kinetic energies, and the second due to and the population difference. We identify the first with a heat conduction 
process and the second with a particle diffusion process. 

One sees by inspection that the choice Na = Nb = N* = N/2 and Ta = Tb = T* represents a symmetric solution 
for all values of the control parameters. The granular temperature of the symmetric state, T*, is given by the nonlinear 
equation: 



1 + na(l 



N* 



(22) 



In order to ascertain the stability of such a symmetric solution we assume Ta = T* + 5Ta, Tb = T* + 5Tb and 
Na — N* + SNa {5Nb — —SNa) and expand the equations to linear order about the fixed point T*,N* with the 
result: 



SNA = -l e -T*/T* 



N*T, 

26Na + 7 ^(6Ta-6Tb) 



(T*y 



(23) 



ST A 



a -T s /T* 



T 

T* 



(ST a -ST b )-(3joj* + 



2 \ 2 
— )ST A 



-T a /T* 



(T* + 2T S ) + 7 w*T* 



SNa 



(24) 



Sf B = -e- T °l T * 

Ta 



2 + ^ 
T* 



(ST A - ST B ) - 37w 



+ - ) 8T B 



2 
~~N* 



_L -TJT' 



(T* + 2T S )+ juj*T* 



SN A 
(25) 

The eigenvalues of the associated 3x3 matrix of coefficients give the three relaxation modes. By adding and 
subtracting the last two equations the linear system factorizes into a decoupled equation for the average value of the 
two temperatures ((STa+STb) /2) with eigenvalue A3 = —(370;* + ^-) and a system of rank 2 involving the temperature 
difference and the occupation number. The more negative of the remaining eigenvalues (we call it A2) is of the same 
order as A3, while the second, say Ai, is smaller in absolute value and vanishes at the "special" temperature, T cr . 
This is obtained by solving the transcendental equation: 



Ts7W . , li. _ I ) - It = 2e~ T °/ T " 

T cr 2 / Tb 



(26) 



Fig. shows that above the temperature T cr , Ai is negative, thus making the symmetric solution stable, whereas 
for temperatures below it Ai is positive. In such a case the symmetric fixed point is unstable and the solution flows 
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away, to an asymmetric fixed point. Hereafter, we shall use the word "critical temperature" or "critical line" instead 
of "special" [Icj. 

Let us also notice that when the typical collision time is much shorter than the characteristic times r c and T b the 
temperature T cr — > 2T s /3, which represents an upper bound for such a quantity. 

At this stage few comments are in order. The transition results from the competition between two effects: the 
diffusion due to the "thermal" agitation of the particles, i.e. the tendency to fill all the available space, and the 
dissipation of energy during collisions which favors clusterization. 

When the external drive is sufficiently weak there appears a second fixed point. This is found by imposing the 
detailed balance: 



N A ~ e - T °/ T B 

and the energy balance 



(27) 



■- (N A T A e- T ° /TA - N B T B e- T ° /TB ) - 2 1 lu a N a T a + -N A (T b — T A ) = (28) 



-(N A T A e~ T °/ TA - N B T B e- T ^ TB )~ 2 1L u B N B T B + -N B (T b — T B ) = 
A good approximation to the highest eigenvalue in such a case reads: 



(29) 



Ai = ^e- T */ T V 7 (T S /T» 



3/2)-f/r b -2e- T */ T *). 



One observes two regions where the solutions are apparently metastable. The first is the low temperature region, 
where a detailed inspection shows that the symmetric state, in reality, is very weakly unstable. From fig. ^we see 
that at low temperatures (T* — * 0) the positive eigenvalue Ai vanishes exponentially as e~ Ts ^ T , thus indicating that 
the symmetric solution might appear stable if the observation time is finite. In fact, the inter-well diffusion for T b — > 
is almost completely suppressed by the Arrhenius factor. 

The second region where the system displays genuine (meta)stability occurs just above T cr . A branch with N A =/= N B 
and T A ^ T B is observed by integrating nmerically the system (12 1H by means of an Euler scheme. These asymmetric 
stationary states are obtained by preparing the system in a subcritical configuration (T ast < T cr ), and successively 
increasing the temperature above T cr . The width of the region is about a ten per cent of T cr . The observed hysteresis 
agrees qualitatively with the results relative to the model of ref . • 

The numerical solutions of the coupled equations is displayed in figs. and We observe that a perturbation 
about the symmetric solution is re- adsorbed for T ast > T cr , whereas for temperatures below T cr the perturbation 
grows initially at an exponential rate, before saturating about a finite value. Let us notice that at low temperatures, 
due to the presence of the Arrhenius factor, the saturation process occurs extremely slowly, a phenomenon which is 
not related to the slowing down which occurs only in the vicinity of T cr . 
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FIG. 2: Evolution of the populations in the two compartments versus time, for three different choices of the heat bath 
temperature, T b = 0.2, 0.3, 1. The remaining parameters are: N = 1000, T s = 1, V = 100, a = 1, Th = 1, T s = 0.5, a = 0.7. 



In fig. 21 we display the line of critical points T cr as a function of the inelasticity (1 — a) for two different values of 
the total number of particles. Above the line the solutions are symmetric (N A — N B ), below are asymmetric. Let us 
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FIG. 3: Granular temperatures in the two compartments corresponding to the evolution of [3 

observe that for a — 1 one finds T cr = 0, because there is no clustering instability in systems of elastic particles. The 
smaller the value of a the higher the value of the critical temperature. We also notice that the slope of the critical 
line is non negative. When the total number of particles decreases, also the critical temperature decreases, because 
of the lower collision rate. 

Figure El illustrates the behavior of the order parameter e = \Na — Nb\/N versus the heat bath temperature for 
two different values of the total number of particles. 

Let us notice that when particles are added to the system the critical point shifts up- wards. This explains why the 
steady number of particles in the less populated compartment decreases when more particles are added 




l-ct 

FIG. 4: Critical line for two different choices of the total number of particles 1000 (below) and 10000 (above) according to 
the mean field theory. Vertical axis temperature, horizontal (1 — a). The line represents the granular critical temperature as 
the inelasticity (1 — a) varies from to 1. Below the critical line the symmetric solution is unstable, above it is stable. The 
remaining parameters are the same as in fig|5| 




T, 



FIG. 5: Asymmetry parameter vs. T;, plotted for two different choices of the total number of particles: 1000 and 1200. The 
remaining parameters are the same as in fig|5| Notice the difference between the two curves and the presence of the hysteresis. 

Finally, we comment that, within our perspective, the approach of Lohse et al. |12| and of Lipowski and Droz [14| is 
equivalent to an adiabatic approximation for the temperature variables. In their cases, the two granular temperatures 
are assumed to be slaved by the occupation number variables. In other words, one postulates that the temperatures 
depend on the instantaneous values of the occupation numbers. We implemented this idea within our approach, 
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but we found a large discrepancy with the previous approximation. In reality, the temperature and the occupation 
variable vary on the same time scale and no slaving principle seems to occur. 



In deriving the evolution equations of the previous section we have implicitly stipulated that the values of the 
occupation number and of the kinetic energy per particle in each box have narrow distributions around their mean 
values and that these quantities are sufficient to characterize the state of the system. Do statistical fluctuations play 
any role or even modify the picture presented above? In the present section, which is not intended to be exhaustive, 
we solve the model of section I by means of the DSMC method. The motivation is twofold: we want to validate the 
picture and the approximations introduced above and unveil some phenomena which are not accounted for by a mean 
field description. 

We simulated two ensembles of Na and Nb particles, respectively, subject to a Gaussian forcing, viscous friction 
and inelastic collisions. In addition, the particles of high energy can change compartment with probability per unit 
time r" 1 . 

The scheme consists of the following ingredients: 

• time is discretized, i.e. t — n ■ dt 

• update all the velocities to simulate the random forcing and the viscous damping: 



where W(t) is a normally distributed deviate with zero mean and unit variance. 

• At every time step and for both compartments, a collision step is performed, i.e. an adequate number of 
randomly chosen pairs of velocities are updated with the collision rule : the pairs are chosen with a probability 
proportional to their relative velocity and normalized in order to have a mean collision frequency as calculated 
in eq. l(T3fl. 

• Place in the other compartment with probability dt/r s particles with kinetic energy greater than T s . 

• Change the time counter n and restart. 

To summarize, at every step each particle experiences a gaussian kick and receives energy from the bath, but 
dissipates a fraction of its kinetic energy by collision and by damping. With respect to previous DSMC simulations 
of granular gases particles can migrate to another compartment, whenever their energy exceeds a fixed threshold. No 
packing effects are included. 

The parameters chosen in the simulation are T 8 = 1, V = 100, a = 1, Tb = 1, t s = 0.5, a — 0.7, whereas the total 
number of particles N and Tj, have been varied. 

First, we analyze the evolution of the order parameter \Na — Nb\/N, when the system is prepared in a symmetric 
configuration. The order parameter displays a behavior similar to that obtained by means of the mean field theory. 
In addition, it displays fluctuations around its asymptotic value. In fig. an average over 200 realizations of the 
evolution of the order parameter \Na ~ Nb\/N is shown. The numerical simulations are in agreement with the 
results of the mean field theory. For > 0.45 the system remains substantially homogeneous. For < 0.45 the 
homogeneous state becomes unstable and the stable configuration that is reached is strongly asymmetric. 

Interestingly, we observe a much slower growth (see inset) near the critical temperature, T cr <~ 0.45. The same kind 
of phenomenon appears at low temperatures (T < T s ), because the transitions from one compartment to the other 
represent very rare events. 

In figure [7| the temperature evolution Ta(£) and Tg(t) (after averaging over 200 realizations and relabeling the 
compartments in such a way that A is always the most populated) are displayed. The plateau values of Ta and Tg 
verify (with small deviations, not larger than 10%) the equation l|27|) . 

The velocity fluctuation can be appreciated by considering the velocity PDF's of the two compartments fig. 
Appreciable deviations from a Gaussian have not been detected observed in all cases considered (neither above, 
neither below, nor in the proximity of the critical temperature). 



MONTE CARLO SIMULATIONS 
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FIG. 6: Evolution of the order parameter versus time. N = 1000, but Tt is varying, while the other parameters are: T s — 1, 
V = 100, a = 1, n = 1, t s = 0.5, a = 0.7. 
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FIG. 7: Temperature evolution for various values of the heat bath temperature The other parameters are the same as in fig. |S| 



The order parameter, instead, shows less trivial fluctuations. A typical trajectory is shown in fig. [21 The associated 
distribution is shown in fig. 1101 We observe that at high temperature the population distribution P(Na) is well 
approximated by a gaussian distribution, but displays two symmetric peaks about the central value Na = 500 when 
the system approaches the critical temperature. Finally, at low temperature we find two very narrow peaks, well 
separated and centered around Na = 950 and Nb — 50. 

Let us consider more closely the behavior of the order parameter near the critical point. According to mean 
field description of section II, at T cr the less negative eigenvalue of the matrix of evolution vanishes. The picture 
looks somehow different when microscopic fluctuations are accounted for. One observes a remarkable behavior in 
the order parameter evolution, namely in the earlier regime it grows in a power law fashion as (Na — Nb) 2 ~ t, 
and not exponentially. To explain this phenomenology we remark that the evolution of the order parameter near 
T cr can be assimilated to that of a particle undergoing a random walk. We assume that the early evolution of 
AN(t) — NA^t) — N/2 can be effectively described by the following equation: 



= j(T b )AN(t) + A(T b ) v (31) 
with 7 lesser than zero above the critical temperature and greater than zero under the critical temperature and r\ 
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FIG. 8: Velocity distribution functions in each compartment for two different values of Tb and TV = 1000. The other parameters 
are the same as in fig. HJ The lines represent Maxwellian velocity distributions, whose variances have been determined from 
fits of the numerical data. 




FIG. 9: Fluctuations of the occupation number NA(t) vs time, when the system is in its statistically stationary state for three 
different choices of Tb- The other parameters are the same as in fig. |S| 



a white gaussian noise. In this case 



C(t, t) = (AN(t + r)AJV(i)) = (ATV(O)ATV(O)) 



27 



, 7 (2t+r) 



2 7 ' 



,7T 



(32) 



which means that if 7 < (asymptotically in a time greater than I/I7I) the equal time correlation function becomes: 



(AN(t)AN(t)) - -f- (33) 



or equivalently A = ^—2~/(AN 2 ), where A7V 2 = lim^oo C(0, t). On the other hand, when 7 = (i.e. 
correspondence of the critical point), the equal time correlation function displays a diffusive behavior: 



(AN{t)AN(t)) -> AH (34) 

that means that the diffusion coefficient for the variable AN is given by D — A 2 / 2. 

Eventually in the late regime the growth saturates due to the onsert of non linear effects. 

In order to verify the plausibility of eq. (|31|l we determine A and 7 by extracting them from numerical measures 
of lim^oo C(r, t) and imposing the asymptotic behavior — ^e 1T . In fig. El we display the values of the constants 
obtained from simulations above the critical temperature. We also added to the plot the value of the diffusion 
coefficient, D = A 2 /2, as measured at the critical temperature. 
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FIG. 10: Probability distribution functions of the occupation numbers in a single well (top and middle panel, respectively) and 
in the two wells (bottom panel) , referred to a system of N = 1000 particles. The other parameters are the same as in fig. 
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FIG. 11: Study of A and 7 versus T b . 



Interestingly, when the characteristic time r s decreases (namely r s = 0.1) a new phenomenon appears: the temper- 
ature of the less populated compartment exhibits an initial regime, during which it increases, due to the decrement 
of the number of collisions the particles experience, followed by a second regime where the temperature instead de- 
creases and eventually reaches a temperature lower than the temperature of the more populated container. This 
kind of anomaly is due to the fact that the fastest and more energetic particles are removed causing a large negative 
fluctuation of the associated average kinetic energy. The positive energy flux provided by the heat bath is not suffi- 
cient to compensate the energy loss due to the removal of the fast particles. This phenomenon is not observed in the 
mean field model, because it results from the large non-gaussian fluctuations of the velocity PDF. We checked the 
fourth cumulant and observed that whereas the cumulant relative to the populated compartment corresponds to a 
gaussian, the cumulant of the small population strongly fluctuates. Experimentally it might be hard to observe such 
a phenomenon, which is probably an artifact of the model, because the limit of small r s is not very realistic. 

Finally, our schematization of the compartmentalized granular gas recalls the Gibbs ensemble method of equilibrium 
statistical mechanics, whereby it is possible to study first order phase coexistence without interfaces |2l|. 



CONCLUSIONS 



To summarize, we introduced a model for a compartmentalized granular gas which allows to bridge between the 
microscopic level to the hydrodynamic level. In the first part of the paper we have derived a Fokker-Planck-Boltzmann 
description starting from the stochastic evolution of the particles coordinates. Next, by employing a gaussian ansatz 
for the velocity distribution function, we have obtained a closed set of equations for the slowly varying fields, namely 
the granular temperatures and occupation numbers of each compartment. Let us comment that with respect to 
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FIG. 12: Populations and granular temperatures of the two compartments when r s is small. Notice that the temperature of 
the less populated compartment initially grows, and successively decreases below the temperature of the other compartment. 

existing flux models, our approach treats analytically the granular temperature on equal footing as the occupation 
variables. The solution of the resulting equations shows the existence of two different regimes: at strong shakings 
the populations in the two compartments are symmetric and for weak shakings the two populations are asymmetric 
together with their granular temperatures. A critical point separates these two behaviors. The dynamics has been 
characterized in the various regimes. In the last part of the paper we have solved the full model by means of DSMC. 
The presence of stochastic fluctuations leads to properties not observed in a mean field description. These are the 
presence of critical fluctuations, of anomalies in the dynamics of the populations and in population fluctuations. 
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